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Raman spectroscopy is a fast, non-destructive means to characterize graphene samples. In par- 
ticular, the Raman spectra are strongly affected by doping. While the change in position and width 
of the G peak can be explained by the non-adiabatic Kohn anomaly at F, the significant doping 
dependence of the 2D peak intensity has not been explained yet. Here we show that this is due 
to a combination of electron-phonon and electron-electron scattering. Under full resonance, the 
photogenerated electron-hole pairs can scatter not just with phonons, but also with doping-induced 
electrons or holes, and this changes the intensity. We explain the doping dependence and show how it 
can be used to determine the corresponding electron-phonon coupling. This is higher than predicted 
by density-functional theory, as a consequence of renormalization by Coulomb interactions. 
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INTRODUCTION 



Graphene is the latest carbon allotrope to be discov- 
ered, and it is now at the center of a significant re- 
search efi^orlji'^ii^i^i^i^. Near-ballistic transport at room 
temperature and high mobility^'^'i'^'Siii^ make it a po- 
tential material for nanoelectronic a^^i^^i^^i^"^ , especially 
for high frequency applications^^. Furthermore, its 
transparency and mechanical properties are ideal for 
micro and nanomechanical systems, thin-film transis- 
tors and transparent and conductive composites and 
electrodesi^'i^'ia. 

Graphene layers can be readily identified in 
terms of number and orientation by inelastic 
and elastic light scattering, such as Raman^fi 
and Rayleigh spectroscopies^L^. Raman spec- 
troscopy also allows monitoring of doping, de- 
fects, strain, disorder, chemical modifications and 

^jg^ j20. 23. 24. 25. 26. 27. 28. 29. 30. 31. 32. 33. 34. 35. 36. 37 Indeed 

Raman spectroscopy is a fast and non-destructive char- 
acterization method for carbons^. They show common 
features in the 800-2000 cm~^ region: the G and D 
peaks, around 1580 and 1350 cm~^, respectively. The 
G peak corresponds to the E2g phonon at the Brillouin 
zone center (F point). The D peak is due to the 
breathing modes of six-atom rings and requires a defect 
for its activation22i22id£. It comes from TO phonons 
around the K point of the Brillouin zone^i^, is active 
by double resonance (DR)^, and is strongly dispersive 
with excitation energy due to a Kohn Anomaly at K— . 
The activation process for the D peak is inter- valley, and 
is shown schematically in Fig. Wl^d): i) a laser induced 
excitation of an electron/hole pair; ii) electron-phonon 
scattering with an exchanged momentum q ~ K; iii) 
defect scattering; iv) electron-hole recombination. DR 
can also happen as intra-valley process, i. e. connecting 
two points belonging to the same cone around K (or 
K'), as shown in Fig. W[i>)- This gives the so-called 
D'peak, which is at'^ 1620 cm~^ in defected graphite 



G peak 



D'peak 



2D 'peak 






D peak 



2D peak 







Figure 1: (Color Online) Role of the electron dispersion (Dirac 
cones, e = ±i)_f|p|, shown by solid black lines) in Raman scat- 
tering: (a) intravalley one-phonon G peak, (b) defect-assisted 
intravalley one-phonon D' peak, (c) intravalley two-phonon 
2D' peak, (d) defect-assisted intervalley one-phonon D peak, 
(e) intervalley two-phonon 2D peak. Vertical solid arrows rep- 
resent interband transitions accompanied by photon absorp- 
tion (blue lines) or emission (red lines) (the photon wavevec- 
tor is neglected). Dashed arrows represent phonon emission. 
Horizontal dotted arrows represent defect scattering. 



measured at 514nm. 

The 2D peak is the second order of the D peak. This 
is a single peak in single layer graphene (SLG), whereas 
it splits in four in bilayer graphene (BLG) , refiecting the 
evolution of the band structure^S. The 2D' peak is the 
second order of the D' peak. Since both 2D and 2D' origi- 
nate from a process where momentum conservation is sat- 



isfied by two phonons with opposite wavevectors (q and 
— q), they do not require the presence of defects for their 
activation, and are thus always present. Indeed, high 
quahty graphene shows the G, 2D and 2D' peaks, but not 
D and D— . Also, under the assumption of electron-hole 
symmetry, the two-phonon peaks are fully resonan t "^^i^^ . 
This means that energy and momentum conservation are 
satisfied in all elementary steps of the Raman process, as 
shown schematically in Fig. [l][c,e). Then, all interme- 
diate electronic states are real. As a consequence, two- 
phonon Raman spectroscopy is sensitive to the dynamics 
of the photo-excited electron-hole pair, in particular, to 
the scattering processes it can undergo. This is of crucial 
importance for the present work. 

The effects of doping on the graphene G-peak position 
[Pos(G)] and Full Width at Half Maximum [FWHM(G)] 
were reported in Refs. [Ulsillsilsi. Pos(G) increases and 
FWHM(G) decreases for both electron and hole doping. 
The G peak stiffening is due to the non-adiabatic removal 
of the Kohn-anomaly at T^iH. The FWHM(G) sharpen- 
ing is due to Pauli blocking of phonon decay into electron- 
hole pairs, when the electron-hole gap is higher than the 
phonon energy^^*^, and saturates for a Fermi shift big- 
ger than half phonon energy2L2^iii. A similar behavior 
is observed for the LO-G~ peak in metallic nanotubea^, 
for the same reasons. In the case of BLG, the differ- 
ent band structure re-normalizes the phonon response to 
doping differently from SLG^'^'^. Also in this case the 
Raman G peak stiffens and sharpens for both electron 
and hole doping, as a result of the non-adiabatic Kohn 
anomaly at F— . However, since BLG has two conduction 
and valence subbands, with splitting dependent on the in- 
terlayer coupling, this changes the slope in the variation 
of Pos(G) with doping, allowing a direct measurement of 
the interlayer coupling strengtb^i^. 

Another significant result is that in SLG the ratio of 
the heights of the 2D and G peaks, I(2D)/I(G), and their 
areas, A(2D)/A(G), is maximum for zero doping^SiiSii^^ 
and decreases for increasing doping. On the other hand, 
this shows little dependence on doping for BLG^>2^. 
Fig. [2l pl ots the combined data for SLG and BLG from 
RefsHiflSllIsi 48 49 50. Note that Refs. reported 
height ratios, while here, as discussed later, we ana- 
lyze the area ratio A(2D) /A(G), which encompasses both 
trends of I(2D)/I(G) and FWHM(2D)/FWHM(G). 

Due to residual disorder, the energy of the Dirac point 
can fiuctuate across the sample on a scale smaller than 
the laser spot, which leads to spatial inhomogeneity of the 
doping level2^>^. We attribute the difference in the be- 
havior of the two SLG curves in Fig. 2 to a different degree 
of residual charge inhomogeneity in the polymeric elec- 
trolyte experiments of Refs. Isiljil . On the other hand, 
the use of this electrolyte enabled probing a very large 
doping range, because the nanometer-thick Debye layer 
gives a much higher gate capacitance compared to the 
usual 300nm Si02 back gate2^>^>^. Note as well that 
A(2D)/A(G) for the most intrinsic samples measured 
to date is ~12-1722iiSii2ii^, much higher than the zero 
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Figure 2: Experimental A(2D)/A(G), measured for 514. 5nm 
excitation, as a function of Ef for slg20,31 ,32,48,49 
BLO^. The BLG data (solid squares) are divided by 10, 
to make comparison easier. Note that the doping dependent 
SLG data are a combination of two experiments on two differ- 
ent samples, from Refls^ (half-filled circles) and Reflsil (Open 
circles), and a data-point representative of intrinsic graphene 
from Refs^MMM (solid star) 
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gating values in Refs. l31l)32l . as shown in Fig. 2. This 
points again to sources of disorder in the gated samples 
of Refs. I3lll32l . while the absence of a significant D peak 
excludes large amounts of structural defects. 

Here, we show that the 2D intensity doping depen- 
dence results from its sensitivity to the scattering of the 
photoexcited electron and hole. Assuming the dominant 
sources of scattering to be phonon emission and electron- 
electron collisions, we note that, while the former is not 
sensitive to doping, the latter is. Then, the 2D doping 
dependence can be used to estimate the corresponding 
electron-phonon coupling (EPC). 



II. DOPING DEPENDENCE OF TWO PHONON 
RAMAN INTENSITY 

A. Theoretical Dependence 

Raman scattering^ is an electron-mediated process 
where electromagnetic radiation exchanges vibrational 
quanta (phonons) with a crystal. A complete description 
requires the detailed knowledge of (i) electronic structure. 
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(ii) phonon dispersions, (iii) mutual interactions between 
electrons and phonons (i.e. electron-electron, electron- 
phonon and phonon-phonon scattering). 

The Raman spectrum of graphene consists of a set of 
distinct peaks. Each characterized by its position width, 
height, and area. The frequency- integrated area under 
each peak represents the probability of the whole process. 
It is more robust with respect to various perturbations 
of the phonon states than width and height. Indeed, 
for an ideal case of dispersionless undamped phonons 
with frequency Wph the shape of the n-phonon peak is a 
Dirac 6 distribution oc 5{oj — nojph), with zero width, in- 
finite height, but well-defined area. If the phonons decay 
(e. g, into other phonons, due to anharmonicity, or into 
electron-hole pairs, due to electron-phonon coupling), the 
S lineshape broadens into a Lorentzian, but the area is 
preserved, as the total number of phonon states cannot be 
changed by such perturbations. If phonons have a weak 
dispersion, states with different momenta contribute at 
slightly different frequencies. This may result in an over- 
all shift and a non-trivial peak shape, but frequency inte- 
gration across the peak means counting all phonon states, 
as in the dispersionless case. Thus, the peak area is pre- 
served, as long as the Raman matrix element itself is 
not changed significantly by the perturbation. The lat- 
ter holds when the perturbation (phonon broadening or 
dispersion) is smaller than the typical energy scale de- 
termining the matrix element. Converting this into a 
time scale using the uncertainty principle we have that, 
if the Raman process is faster than the phonon decay, 
the total number of photons emitted within a given peak 
(i. e., integrated over frequency across the peak), is not 
affected by phonon decay, although their spectral dis- 
tribution can be. Although the graphene phonons giv- 
ing rise to the D and D' peak are dispersive due to the 
Kohn Anomalies at K and F— , their relative change 
with respect to the average phonon energy is at most 
a few %, thus we are in the weakly dispersive case dis- 
cussed above. The phonon decay in graphene is in the 
picosecond timescale, while the Raman process is faster, 
in the femtosecond timescal o^^i^^i^^ . Then, we will an- 
alyze the area ratio, A(2D)/A(G), which encompasses 
both variations in height ratio, I(2D)/I(G), and width: 
FWHM(2D)/FWHM(G). 

We first consider the G peak. For the one-phonon pro- 
cess, allowed by momentum conservation, which gives rise 
to the G peak, the picture is entirely different from the 
two-phonon case. As shown in Fig. la, the process respon- 
sible for the G peak is determined by virtual electron- 
hole pairs with energy El/2, where El is the laser exci- 
tation energy (for a typical visible Raman measurement 
El/2 ~ leV). If the Fermi energy, Ep, stays below El/2, 
as in Refs. [3lll32 . these electronic states are not strongly 
affected. Only the final phonon state is infiuenced by 
doping, which manifests itself in a change of Pos(G) and 
FWHM(G)^>2ii^i2^. However, the area of the peak is de- 
termined by the total spectral weight of the phonon state, 
which is preserved. Thus, we do not expect any signifi- 



cant dependence of A(G) on doping, as long as the dop- 
ing is not too strong, so that \Ep\ 1 eV. We can then 
take the measured doping dependence of A(2D)/A(G) as 
representative of the A(2D) trend. Note that A(G) can 
change as a function of other external parameters, such 
as the Raman excitation energy^Si^'^i^'^. However, 
for fixed excitation, such as in the experiments discussed 
here, the above argument holds. 

In RefSi the following expressions for the 2D and 2D' 
areas were obtained: 




where e is the electron charge, c is the speed of light, 
e^/c « 1/137 is the fine structure constant, and vf 
is the electron velocity (its experimental value is vp 
10^ m/s « 6.6 eV ■ }M>M^). 2j is the scattering rate of 
the photoexcited electron and hole. Note that we define 
7 as the imaginary part of the energy, so it determines 
the decay of the amplitude, while the decay of the prob- 
ability is determined by 2j. This includes all sources of 
inelastic scattering. Assuming the two main mechanisms 
for electron scattering to be the emission of phonons and 
electron-electron collisions, we write: 

7 = 7o-ph + 7oe, 7e-ph = 7r + 7k- (2) 

Here we include the phonons near T and K, responsible 
for D and D'. The corresponding emission rates, 27r and 
2'yK, enter the numerators in Eqs. I|lap . (|lb|) . 

Two points regarding Eqs. (fTa|) . (HbJ should be em- 
phasized. First, the scattering rates depend on the elec- 
tron energy, e, which is defined by half the laser energy, 
e w El/2 [see Eq. I|12p in the next section]. Second, if 
impurity scattering is significant compared to other scat- 
tering mechanisms, the corresponding elastic scattering 
rate cannot be simply included in 7 and Eqs. lfTa|) . l|lbp . 
The whole Raman intensity calculation should be done 
differently. Eqs. (fTa|) . thus neglect impurity scatter- 
ing. For short-range impurities this assumption is justi- 
fied by the absence of a large D peak in the spectra of 
Refs. fsills^ . Long-range disorder is efficiently screened 
(even though the vanishing density of states at the Dirac 
point requires the screening to be nonlinea r^^i^^i^^i^^ ); 
it is precisely this screening that gives rise to the inho- 
mogeneous concentration of electrons/holes and spatial 
fiuctuations of the Dirac point energy. 

In principle, there are no reasons for a strong depen- 
dence of 7c-ph on carrier density. However, 7oe does ex- 
hibit such a dependence. Indeed, in undoped graphene at 
low temperatures, the photoexcited electron finds itself 
in a state with some momentum, p, measured from the 
Dirac point, in the empty conduction band. To scatter 
into a state with a different momentum p', it has to give 
away some energy and momentum to another electron in 
the full valence band. This second electron would have 
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to be promoted to the conduction band (as there are no 
available empty states in the valence band) into a state 
with momentum Pe, leaving a hole in the valence band 
with ph- Momentum and energy conservation require: 

P = p' + Pe + P/i, (3a) 
e(p) = e(p') + e(Pe) + e(Ph), (3b) 

where e(p) is the quasiparticle dispersion, assumed the 
same for electrons and holes. For Dirac particles, e(p) = 
wfIpI, the only possibility to satisfy both conservation 
laws is to have all four momenta parallel. If the spectrum 
is convex, (Pe{p)/dp^ > 0, the two equations can be sat- 
isfied by a set of momenta with non-zero measure, i. e. 
the phase space is finite. If it is concave, (Pe{p)/dp^ < 0, 
they are incompatible. In SLG the spectrum is Dirac to 
a first approximation, resulting in an uncertaintj*^. This 
can be resolved by taking into account corrections from 
electron-electron interactions, which make the spectrum 
concave , ^^1^'' and the interband process forbidden. 

As new carriers are added to the system, intraband 
electron-electron collisions become allowed. The momen- 
tum and energy conservation become: 



P + Pe = P 

e(p) + e(Pe) 



^(p')+<Pe), 



(4) 
(5) 



which can be satisfied for any quasiparticle dispersion. 
These collisions give a contribution to 7ee which increases 
with carrier concentration. As a consequence, the to- 
tal 7 in Eq.la increases, leading to an overall decrease of 
A(2D), consistent with the experimental trend in Fig.[2l 

The above arguments essentially use the non-convexity 
of the electronic spectrum in the conduction band, and 
thus apply to SLG only. In BLG, the spectrum is 
parabolic near the Dirac point, so that cPe/dp^ > 0, and 
the phase-space restrictions are absent. Thus, electron- 
electron collisions are allowed even at zero doping, and 
the collision rate has a much weaker dependence on Ep, 
which, in first approximation, can be neglected. Thus, 
A(2D) is expected to have a weak dependence on Ep, as 
seen in Fig. [21 where the experimental A(2D)/A(G) for 
BLG shows a negligible variation with doping^. 

To quantify the doping effects on the SLG A(2D), 
we first calculate the electron-electron scattering rate, 
27oo, in the random-phase approximation, analogously to 
Refs. IssIIgqI . 7ee is given by the imaginary part of the on- 
shell electronic self-energy, ImHeeip, e) for e — *■ vpp — O^, 
with e and p counted from the Dirac poinli^. Here we 
consider the limiting case, when the energy of the pho- 
toexcited electron (e = El/2) far exceeds Ep. The car- 
rier concentration is n = Ep/{Trvp'^). In this case, the 
collisions are dominated by small momentum transfers, 
|P — P'l ~ \Ep\/vp, so 7oo does not depend on e and 
is proportional to \Ep\, the proportionality coefficient 
depending only on the dimensionless Coulomb coupling 
constant = e^/{evp) (e being the dielectric constant): 



7e 



= \Ep 



\svp 



+ OiEl/e), 



(6) 




Figure 3: Numerical values of /(rs), from Eq.® 
where the function / is given by: 

firs) j 





2/(l+cos<p) 



dx sin ip Ri 







2/(1— cos <p) 



[2{x/rs + 4)x sin ip]^ + Rj 



dx x'^ sin ip i?2 



[2(x/r, + 4)x sin <f - R^f Rl{x, >p) 



2/(l+cosi^) 

and Ri, i?2, ^3 are: 

Ri{x, ip) = a+6_|_ — a_fe_ — x'^ In 
R2{x, (p) — a+6+ — x'^ In 



(7) 

, (8a) 
(8b) 

Rsix, ip) = a- J x^ — a}_ — x^ arccos — , (8c) 



2 «+ + fe+ 
a- + 6_ 



2 ,„ a+ + 



a± = 2 ± a; cos v?, 6± = y a± ~ ■ 

Fig. [3] plots f{rs), calculated numerically. 
Thus, we expect A(2D) to change with Ep as: 



(8d) 



A{2D) 



C 



[7e„p„ + \Ep\f{eyevp)]^ 



(9) 



with C a constant. Note that a variation of the dielectric 
constant e will affect A(2D). Given the negligible depen- 
dence of A(G) on doping, Eq. ([9]) can be rewritten as 



AiG) 



V M'2D) 
where C is another constant. 



= C"[7e-ph + \Ep\f{eyevp)], (10) 
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Figure 4: Fit of the experimental dependence \/ A(G) /A{2D) 
from Ref. [3l| (open circles) and Ref. [H (half-filled circles) 
using Eq. I|lip (dashed and solid lines, respectively). 



B. Fit to Experiments 



Fig. H plots .JA{G)/A{2D) as a function of Ep. This 
dependence, according to Eq. (flOl) . should correspond 
to two symmetric straight lines joining at Ep = 0. As 
noted in Sec. I, close to Ep = the data from the two 
polymer electrolyte gating experiments do not converge 
to the same va lue. However, for both a linear rise of 
y/ A{G)/A{2D) is seen at higher energies. Also, while 
the data represented by open circles in Fig. 4 are almost 
symmetric, a significant asymmetry is seen for electron 
doping in the set represented by the half-filled circles, 
while the two sets are in good agreement for hole doping. 

A(2D)/A(G) for intrinsic samples measured at 514.5 
nm excitation, the same used in Refs lsUfsl . is in the 
range 12-l?22ii2iii£, represented by the star in Fig. 4 at 
14.5. This is in good agreement with the ratio mea- 
sured for carbon whiskers^. These show a 2D peak 
very similar to graphene, being composed of mis-oriented 
graphene layer a^°i^° . However, their Raman spectra are 
much less susceptible to charged impurities or surface 
d oping, being bulk materials^. This corresponds to 
y/A{G)/A{2D) - 0.24 - 0.29, which we use to elimi- 
n ate the effect o f doping inhomogeneity, by constraining 
,yA{G)/A{2D) - 0.26 at zero doping. We also need 
to consider the dielectric constant of the polymer elec- 
trolyte^i, e = 5, giving f{e'^/evp) « 0.06. Thus, we fit 



the data with a one-parameter expression 
0.26 



I AjG) 

A{2D) 



(7e-p/. + 0.06|eF|). (11) 



We fit separately each branch of the two data-sets, as 
shown by solid and dotted lines in Fig. [H We get je-ph- 
18, 21, 29, 65 meV, with an average je-ph ~ 33 meV. 



III. RAMAN INTENSITIES AND 
ELECTRON-PHONON COUPLING 

A. Theoretical Background and Electron-Phonon 
Coupling Definitions 

Even though graphite and other sp^ -hybridized mate- 
rials have been investigated for more than 50 years^SiZi, 
all the fundamental physical properties needed for the 
interpretation of the Raman spectra have undergone an 
intense debate, which seems to be just beginning to con- 
verge. Interestingly, several features of both phonon 
dispersions and band structure of graphene are deter- 
mined by the EPC. For example, in the Kohn anomalies 
around T or K— the correction to the phonon frequen- 
cies due to EPC results in a linear slope of the opti- 
cal phonon branches as the wave vector approaches T 
or K. The EPC and phonon dispersions calculations of 
Ref. m have been confirmed at the F point by inelastic X- 
ray scattering^^, and by the measurement of FWHM(G) 
in graphite, graphene and nanotubes22ii2^>^i^, once an- 
harmonic effects are taken into accoun t^"'^^'^^ . For the K 
point, the precise slope of the anomaly is debated^^i^liZ^. 
Another EPC effect is the kink in the electron dispersion, 
~200 meV below Ep, seen by angle- resolved photoemis- 
sion spectroscopy (ARPES)^>2^. This is attributed to a 
correction to the electron energy due to EPG^'^'21, al- 
though alternative explanations also exis1j2^. Thus, a cor- 
rect EPC determination is a fundamental step for an ac- 
curate description of the physical properties of graphene, 
and nanotubes, being rolled up graphene sheets. 

To link the 2D intensity to the EPC we first consider 
the rate of phonon emission by the photoexcited elec- 
tron/hole, 27o_ph. This is obtained from the imaginary 
part of the electron self-energy, 7o-ph = IniSc-ph(e)- For 
El/2 > Ep + ujr, as in the case of the Raman measure- 
ments at 2.41 eV excitation of Refsl3lll32. we have^: 



7-ftr 



4 



El 
2 



Ar 

7r = — 



El 

2 



Then, from Eq. (2): 











7c-ph = — 









(12) 



(13) 



The dimensionless coupHng constants Ar, \k correspond 
to phonons close to T and K, respectively, and determine 
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their rate of emission. We define them as: 



f'^ a 



2MuJY,KVp 



(14) 



Here ujk = 1210 cm-^ = 0.150 eVl^ and lot = 
1580 cm-i = 0.196 eV,^ M « 2.00 • 10"^^ g = 2.88 ■ 
10^ (eV • A is the mass of the carbon atom, ^u.c. ~ 
5.24 A is the unit cell area. Ft and Fk have the dimen- 
sionality of a force and are the proportionality coefficients 
between the change in effective hamiltonian and the lat- 
tice displacement along the corresponding phonon mode. 
Strictly speaking, the relevant phonon states are not ex- 
actly at r and K, as shown in Fig. 1. However, the 
corresponding deviation, q ^ El/vf, is small compared 
to the K-K' distance, and is neglected. All observables 
depend on the dimensionless EPCs, Ar and Xk- 

Eq. ([H]) follows the notation of Ref. [13. Since dif- 
ferent EPC definitions are used in the literature, it is 
quite useful to give here matching rules for all of them, 
which will be necessary when comparing the EPC val- 
ues obtained here with previous (and future) reports. 
The EPCs can be conveniently matched by either re- 
lating them to the nearest-neighbor tight-binding model, 
where the constants are expressed in terms of a single pa- 
rameter: dto/da, the derivative of the nearest-neighbor 
electronic matrix element with respect to the interatomic 
distance, or by comparing expressions for various observ- 
ables. For example, doping leads to a G peak shift due 
to EPC. This is expressed in terms of Ep as^^i^i^: 



^'^r = ^ ( \Ef\ 



— m 

4 2Ef + ujr 



(15) 



The corrections to the phonon dispersions as function of 
wavevector q, measured from F or K, are^^i^i^^: 



/ 2 9 2 

n'^ — 



5uj 



r-TO 



Ar cjp 



8 y^vj.q'^ - UJ^ 



(16a) 
(16b) 

(16c) 



Note that the i?2g mode splits into longitudinal (F — LO) 
and transverse (F — TO) at finite q. Note also that 
due to analytical properties of the logarithm and square 
root, Eq. ^ at \Ef\ < wr/2 and Eqs. (fT6a|) - fT6c|) at 
vpq < ujk.t acquire imaginary parts, which correspond 
to the phonon decaying into a continuum of electron-hole 
pairs^. In this case 2 Im Sui gives the FWHM of the cor- 
responding Lorentzian profile. At vpq ^ i^K,r Eqs. I|16ap 
and l|16cp give the profile of the Kohn anomalies. 

In Refs.lSli^Ilillithe EPCs are defined as the matrix 
elements of the Kohn-Sham potential, differentiated with 
respect to the phonon displacements. What enters the 
observables are their squares, averaged over the Fermi 



surface in the limit Ep 0. The matching rule is then: 

(17a) 
(17b) 

In Ref. m the dimensionless coupling constant A is de- 
fined as the proportionality coefficient in Eq. (fT5|) . Thus, 



27r' 



(18) 



Note that the expression linking EPC to FWHM(G) in 
Ref. m underestimates FWHM(G) by a factor 2. 

The dimensionless EPC reported in the ARPES anal- 
ysis of Refs. liollT^lsollsil and in the scanning tunneling 
spectroscopy (STS) experiment of Ref. [H was measured 
from the ratio of the electronic velocities below and above 
the kink in the electron dispersion. This ratio is deter- 
mined by the derivative of the real part of the electronic 
self-energy ReSe-ph(e) due to the EPC. The latter can 
be calculated if one takes the Dirac spectrum for electrons 
and a constant dispersion for phonons. For Ep > 0^: 

Se-p/i(e) = - — (e-tJA-)lnr — ^ 

47r I e — ujii — Ep\ 

Xk , , N, EM\£ + i^K — Ep\ 

- — e + ujk) In — -2 

47r (e -I- ujk) 

Ar / M Em 

47r \e — — Ep\ 

Ar, , ^-, Em\(- + <^t - Ep\ 

- — e + wr) In — . 19 

47r (e + ojr) 

Here Em is the ultraviolet cutoff, of the order of the 
electronic bandwidth. We then get the matching rule: 



^(kink) 



9RcI]e- 



Xk f Ep 
2^ V 



<i = EB 



LOK 



In- 



Em 



UJK 



UK + Ep 



At 
27r 



Ep — ujy 
uJr 



In- 



Em 



UJr + Ep 



(20) 



However, we note that Xk is subject to Coulomb 
renormalizationsS^. This implies that Xk depends on the 
electronic energy scale, such as the electron energy e, the 
Fermi energy Ep, or the temperature T, whichever is 
larger: Xk = Aj^:(max||e|, |, Tj). This dependence is 
shown in Fig. 6 of Ref. |83| in the semi-logarithmic scale. 
In a Raman measurement this scale is given by the en- 
ergy of the photo-excited electron: e « El/2, as long 
as El /2 > \Ep\. Thus, in Eq. ^ Xk = Xk{El/2). 
On the other hand, to estimate the EPC effects on the 
phonon dispersions in the intrinsic graphene, the rele- 
vant electron energy is of the order of the phonon energy. 
Thus, in Eq. (fTBc]) Xk ~ Xk{ujk)- From Fig. 6 of Ref.tsa 
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we estimate that \k{'^k)I^k{ElI'2) ~ 1-5 for e = \ and 
1.2 for e = 5 (taking « 2 eV to represent Raman 
measurements in the visible range). 

The situation with Eq. ((20l) is more compHcated, since 
the cutoff A/ appears expHcitly. The logarithmic term is 
determined by all energy scales from Ej^i down to Ep + 
ujk- Thus, the proper expression is 



^(kink) 



^k{Ef) Ep — luk 
27r ujk 

^ Ar / Ep — loy 
2tt \ ujr 



27r e 



In 



Em 
Ep + ujr 



(21) 



B. Experimental Electron-Phonon Coupling 

From Eq.(12), our overall average je-ph = 33 meV, 
derived from a fit to all the data in Fig. 4, gives: 



Ar + A 



K 



0.13. 



(22) 



On the other hand, the hole doping side of Fig. 4 shows 
two data sets very consistent with each other. We can 
thus get another estimate taken from the average je-ph ~ 
20 meV for just the hole doping side. This would give: 



Ar + A 



K 



0.08. 



(23) 



Based on measurement a^^i^^ and DFT calculations^^, 
the value of Ar can be reliably taken « 0.03. Indeed, 
DFT gives2£ {gl)p = 0.0405eV^ and vp = 5.5eV- A, cor- 
responding, from Eqs. (HI]), l|17ap to Ar ~ 0.028. Even 
though {gi')F and vp are subject to Coulomb renormal- 
ization, Ar = 'iAu.c.{9r) f /vp, which contains their ratio, 
is not.— The experimental Ar extracted from FWHM(G) 
in graphene and graphite^^iii according to Eq. (|16ap and 
from the dependence of Pos(G) on Fermi energy accord- 
ing to Eq. mi), give Ar « 0.034^^ and Ar « 0.027-2^. 

On the other hand, the value of Xk is still 
debated2i>2a>M. The calculated DFT {g^)p = 0.0994eV^ 
together with the DFT vp = 5.5 eV ■ A (both taken from 
Ref. lli) gives Xk = 0.034. However, Ref. HI suggested 
this should be enhanced by Coulomb renormalization by 
up to a factor 3, depending on the background dielec- 
tric constant. In order to compare with our fits, we need 
consider that the corrections to the phonon dispersion are 
determined by electronic states with energies lower than 
those contributing to the Raman signal. As discussed in 
Sec. IIIA, Xk{ujk)/Xk{El/2) w 1.2 for e = 5. Our fit in 
Eq. ((221) corresponds to Xk{El/2) w 0.1, while Eq. ^ 
gives Xk{El/2) « 0.05, resulting in Xk{^^k) ~ 0.12 and 
Xxii^K) ~ 0.06, respectively. These are bigger than DFT 
by a factor of about 3.5 and 1.7, respectively. 



A recent GW calculation gave {D^)p = 193eVVA' 



74 



Combining this with the GW vp = 6.6 eV ■ Aj^i, we get 
Xxii^K) ~ 0.054, a factor ^ 1.6 greater than DFT, in 
good agreement with our fitted average on the hole side. 



Ref. [73 reported inelastic x-ray scattering measure- 
ments of the phonon dispersions near K more detailed 
than those originally done in Ref. [tI, now giving a 
phonon slope at K of 73 meV • A. Using Eq. I|16cp at 
q ^ ujk/vp and taking the experimental value vp = 
6.6 eV ■ A^ (the bare electron velocity, i. e. below the 
phonon kink), we obtain Xxi^jK) ~ 0.044, a factor ^ 1.3 
higher than DFT, again in good agreement with our fit- 
ted average on the hole side. 

Another EEC estimate can be derived from the 2D and 
2D' area ratio. Combining Eqs. (fTa|) . (fTbl) .12.13 we get: 



A{2D) 
A{2D') 



K 



(24) 



For intrinsic SLG and graphite whiskers, the experimen- 
tal A(2D)/A(2D') is - 25 - 302£^^^, which gives 
Xk{El/2) « 0.11 and Ar + Xk{El/2) w 0.13. Since 
in this case s — 1, this results in Xk{u!k) ~ 0.16, a fac- 
tor ~4.5 higher than DFT, in agreement with our upper 
estimate from Eq. (22). 

We finally consider the EPC derived from ARPES and 
STS. For an estimate, we approximate the dependence 
Xk{£) as Hnear in Ine. We take Xk{Em) = {^tI^k)Xt, 
as given by DFT (assumed to be valid at high energies), 
and leave Xk{El/2 ~ 1 eV) as the only free parameter 
determining this Hnear dependence: 



Aif(e) = — Ar- 

UJK 



UJr 

— Xt-Xk{El/2) 

UK 



Taking Ep = 0.4 eV^^^^, Eaj = 
tuting Eq. ^ in Eq. we get: 



HEn/e) 
\n[EM/{EL/2)y 
(25) 

10 eV, and substi- 



;^(ki„k) «o.7Ar + 0.6Aif(i;L/2). 



(26) 



Note that the dependence on the precise value of Em 
is weak: setting Em = 5 eV changes the first coeffi- 
cient to 0.5, and the second (more important as it mul- 
tiplies the larger coupling constant) varies only by 2%. 
The measurements in Refs. %0,;i6,MMM gave A'''"*) w 
0.4, 0.3, 0.26, 0.2, 0.14, respectively. The smallest of 
these values, A^'*^'"''^ ~ 0.14, from Eq. ((26l) corresponds 
to Ar + Xk{El/2) « 0.23, while the highest to Ar + 
Xk{El/'2) ~ 0.66. Even the smallest is almost twice our 
upper bound fit of Eq. ((22l) and would imply an EPC 
renormalization of almost one order of magnitude. Res- 
olution effects could play a role in this overestimatior^. 

Thus, our fits to the doping dependent Raman area 
ratios point to a significant renormahsation, by a factor 
1.7-3.5, of the TO mode close to K, responsible for the 
Raman D and 2D peaks. Our lower bound estimate is 
consistent with recent GW calculations and phonon mea- 
surements, but our upper bound is much lower than the 
smallest estimate derived by ARPES. 
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IV. CONCLUSIONS 

We have shown that the 2D intensity dependence on 
doping can be explained considering the influence of 
electron-electron interactions on the total scattering rate 
of the photogenerated electrons (holes) . We have given 
a simple formula linking 2D peak area to the Fermi level 
shift. Fitting this to the available experimental data we 
got an estimate for the EPC value of the TO phonons 
close to K, responsible for the Raman D and 2D peaks. 
This is larger than that from DFT calculations, due to 



renormalisation by Coulomb interactions. However, our 
fitted EPC is still significantly smaller than those re- 
ported in ARPES or STS experiments. 
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